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A  capillary  discharge  based  pulsed  plasma  thruster  is  currently  under  development  at 
the  Air  Force  Research  Laboratory.  A  torsion  thrust  stand  has  been  developed  to 
simultaneously  measure  the  shot  by  shot  total  impulse  and  propellant  mass  loss  of  the 
thruster,  yielding  a  per  shot  specific  impulse.  Competing  design  considerations  affecting  the 
performance  of  the  thrust  stand  has  lead  to  an  optimized  diagnostic  tool  capable  of 
measuring  total  impulses  up  to  150  mN  s  with  a  resolution  of  0.4  mN  s  and  mass  losses  up  to 
18g  with  a  resolution  of  1  mg. 
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torsion  damping  coefficient  [N  m  s/rad] 

logarithmic  decrement  [n.d.] 

damping  ratio  [n.d.] 

steady  state  calibration  force  [N] 

gravitational  constant  [m/s2] 

mass  moment  of  inertia  [N  m  S2/rad]  [kg  im/rad] 

total  impulse  [N  s] 

specific  impulse  [s] 

torsion  spring  constant  [N  m/rad] 

mass  of  propellant  [kg] 

radial  distance  from  pivot  to  calibration  force  [m] 
radial  distance  from  pivot  to  translation  sensor  [m] 
radial  distance  from  pivot  to  CD  [m] 
rotational  deflection  of  the  thrust  stand  [rad] 

linear  displacement  of  the  thrust  stand,  at  the  displacement  sensor  [m] 
natural  frequency  of  the  thrust  stand  [rad/s] 
damped  frequency  of  the  thrust  stand  [rad/s] 


I.  Introduction 

A  capillary  discharge  based  coaxial,  electrothermal  pulsed  plasma  thruster  (PPT)  is  currently  under  development 
at  the  Air  Force  Research  Laboratory  (AFRL)  as  a  high  efficiency  alternative  to  more  traditional  ablative  PPTs. 
The  major  drawback  of  the  traditional  ablative  PPT  is  the  propulsive  efficiency.  Even  though  specific  impulses  (Isp) 
on  the  order  of  1000  seconds  have  been  achieved  with  ablative  PPTs  utilizing  electromagnetic  acceleration,  the 
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propulsive  efficiencies  of  flight  qualified  thrusters  typically  remains  below  10%. 1  Capillary  discharges  (CD)  are 
relatively  efficient  sources  of  high  density,  high  temperature  plasmas  which  are  being  developed  for  a  number  of 
applications.2,3  4  The  capillary  discharge  being  developed  at  AFRL  consists  of  a  high-density  polyethylene  (HDPE) 
capillary  tube,  a  tungsten  anode  at  the  closed  end  of  the  tube,  and  a  stainless- steel  cathode  at  the  open  end.  The 
anode  is  multi -segmented  with  differentially  applied  potentials  to  provide  ignition.  Capillary  discharges  maintain  a 
resistive  arc  (generated  at  the  anode)  through  the  insulating  capillary  by  the  continual  ablation  of  the  HDPE  wall 
material.  Previous  performance  calculations  have  suggested  that  a  capillary  discharge  based  electrothermal  PPT  can 
achieve  propulsive  efficiencies  of  30-40%  even  without  nozzle  expansion.5  Additional  studies  indicate  that  nozzle 
expansion  of  the  high  pressure,  high  temperature  plasma  generated  by  the  capillary  discharge  will  increase  the 
efficiency  above  50%.6  With  the  potential  for  vast  improvement  in  efficiency,  the  capillary  discharge, 
electrothermal  PPT  may  be  a  strong  candidate  for  satellite  orbital  insertion,  attitude  control,  or  station-keeping 
maneuvers. 

In  the  present  study,  a  diagnostic  tool  has  been  developed  to  investigate  the  applicability  of  a  capillary  discharge 
as  an  electrothermal  PPT.  A  torsion  thrust  stand  similar  to  that  described  by  Ketsdever,  et  al.7  has  been  designed  to 
simultaneously  measure  the  impulse  and  mass  loss  of  the  capillary  discharge.  The  simultaneous  nature  of  the 
developed  technique  will  allow  a  per-pulse  measurement  of  the  discharge’s  specific  impulse  through  the  relation 


(1) 


The  technique  is  especially  useful  for  propulsion  systems  that  utilize  solid  propellant  where  a  direct  measurement  of 
the  propellant  mass  flow  is  extremely  difficult.  The  thruster  will  be  configured  on  the  thrust  stand  such  that  the 
impulse  generated  by  the  discharge  and  the  steady-state  force  generated  by  the  propellant  mass  loss  act  in  the  same 
direction.  The  combined  signal  from  these  effects  can  then  be  decoupled  to  assess  the  ratio  of  the  impulse  to  the 
weight  of  propellant  expended  yielding  the  specific  impulse.  The  optimization  for  the  thrust  stand  will  be  discussed 
in  following  sections. 

A  second  thrust  balance  will  then  be  configured  in  a  manner  similar  to  that  described  by  D’ Souza  and  Ketsdever8 
such  that  the  stand  will  measure  the  time  accurate  shape  and  magnitude  of  the  impulse  delivered  by  the  capillary 
discharge.  Previous  experimental  measurements  have  shown  that  the  typical  capillary  discharge  pulse  lasts  several 
hundred  micro-seconds.5  By  investigating  the  time  accurate  pulse  by  pulse  impulse  delivered  by  the  capillary 
discharge,  discharge  efficiency  may  be  dramatically  improved.  The  time  resolved  impulse  measurement  (TRIM) 
technique  involves  investigating  the  equations  of  motion  for  the  torsion  thrust  stand  and  recreating  the  forcing 
function  from  the  measured  time  resolved  displacement  of  the  stand  as  described  herein. 


II.  Thrust  Stand  Model  and  Analysis 

The  thrust  stand  was  designed  as  a  damped  torsion  force  balance.  The  motion  of  the  thaist  stand  can  be 
modelled  as  an  under-damped  second  order  system  described  by 

70(0  +  C0(O  +  KQ(t)  =  M(t)  (2) 

where  M(t)  is  the  forcing  moment.  In  order  to  model  the  expected  motion  of  a  potential  design,  eqn.  (2)  can  be 
solved  for  0  such  that9 
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where  M0  is  a  steady  state  forcing  moment.  For  time  dependant  forcing  moments,  these  equations  can  be  discretized 
in  time  and  solved  over  individual  time  steps  forced  by  a  constant  forcing  moment, M0  =M(t),t  =  t0  —>t0+At. 

Using  the  small  angle  approximation,  sin 9^6  and  the  geometric  relations  0=S/R,  where  S  is  the  arc  swept  by  a  point 
on  the  thrust  stand,  and  sin 0=XZR,  where  X  is  the  linear  displacement  of  a  point  on  the  thrust  stand,  the  rotational 
equations  can  be  linearized  at  a  radial  distance,  R.  This  yields  a  relation  from  eqn.  (2)  in  terms  of  the  signal  from  a 
displacement  detector  and  the  forcing  function  of 
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As  shown  in  Ref.  7,  the  simultaneous  measurement  of  impulse  and  mass  loss,  where  the  ratio  yields  Isp,  is 
conducted  by  analyzing  two  aspects  of  the  thrust  stand's  overall  motion.  The  maximum  deflection  range  of  the 
thrust  stand  during  a  test  is  linear  with  the  magnitude  of  the  impulse.  The  change  in  the  thrust  stand  position  before 
and  after  the  impulse  applied  is  linear  with  the  propellant  mass  loss.  An  example  of  these  two  measurements  can  be 
seen  in  Figure  1.  The  assumptions  necessary  for  the  linear  relations  to  hold  are  that  the  pulse  duration  be  much  less 
than  the  stand’s  period  and  the  mass  loss  has  a  known  distribution  as  a  function  of  radial  distance  from  the  pivot.  If 
these  conditions  are  satisfied,  two  sets  of  calibration  runs,  known  impulse  versus  deflection  range  and  known  steady 
state  force  versus  change  in  zero,  are  required  to  extract  the  absolute  impulse  magnitude  and  total  mass  loss  from  a 
single  thrust  stand  trace. 
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Figure  1  Simulated  trace  of  expected  impulse  and  mass  loss  for  the  operational  thruster. 

For  the  measurement  of  the  time  dependant  impulse  shape  (profile)  using  the  TRIM  method,  Ref.  8  shows  the 
validity  of  measuring  the  left  hand  side  of  eqn.  (6)  to  reproduce  the  right  hand  side.  In  order  to  accomplish  this,  the 
constants,  K ,  C,  and  /,  must  first  be  determined  for  the  thrust  stand.  To  find  K ,  a  steady  state  force  is  applied  to  the 
thrust  stand.  Given  sufficient  time  to  allow  for  the  damping  of  oscillations,  the  limit  of  eqn.  (6)  as  ►  oc  becomes 


K^FR 
R  c  c 
s 


(7) 


The  remaining  two  coefficients  can  be  found  by  analyzing  the  measured  period  and  the  decay  of  the  deflection 

amplitude  as  a  function  of  time.  First,  a  logarithmic  decrement  is  defined  as 

8  =\n(X./X.  )  (8) 

n  v  i  i+n  v  ' 

where  i  is  a  positive  integer,  and  Xt  and  Xi+n  are  peak  amplitudes,  separated  by  n  periods.  The  damping  ratio  is  then 

given  by 


(n-2n)2+8  2 
n 


(9) 


The  measured  frequency  of  the  stand  can  be  found  from  the  measured  period  and  related  to  the  natural  frequency  by 

(10) 

where  T  is  the  measured  period  of  the  stand,  and  the  mass  moment  of  inertia  and  damping  coefficients  can  be 
calculated  as 


I  =  ~2  (U) 

% 

C  =  2<z4kI  (12) 

Therefore,  the  coefficients  I,  C,  and  K  may  be  determined  from  any  single  test  trace,  or  averaged  over  many  test 
traces,  resulting  from  a  known,  static,  steady-state,  calibration  force.  It  follows  that  if  the  deflection,  X(t %  resulting 
from  a  dynamic  force,  is  measured  on  a  stand  with  no  significant  change  in  I,  C,  or  K  during  that  time,  then  the  time 
dependent  function  F(t)can  be  determined  through  equation  (6). 
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III.  Thrust  Stand  Design 

The  major  components  of  the  thrust  stand  are  the  rotational  arm,  damping  magnets,  flexure  pivots  (torsion 
springs),  displacement  sensor,  and  calibration  system.  The  rotational  arm  contributes  to  the  moment  of  inertia  for 
the  stand  and  provides  a  structure  to  attach  thruster  other  components.  The  damping  magnets  are  used  as  eddy 
current  generators  to  damp  the  stand’s  motion  and  define  the  damping  coefficient  for  the  stand.  The  torsion  springs 
contribute  to  the  spring  constant  by  providing  the  restoring  force.  The  displacement  sensor  measures  the  linear 
motion  of  the  stand  and  defines  the  limits  of  the  stand’s  motion  as  well  as  the  precision  of  the  measurements.  The 
calibration  system  provides  a  known  steady  state  force  and  transient  impulse  to  the  stand  and  defines  the  useful 
range  of  the  stand’s  accuracy.  The  major  components  of  the  stand  can  be  seen  in  Figure  2. 


Figure  2  CAD  drawing  of  first  iteration  thrust  stand  in  the  impulse  profile  measurement  configuration. 

Special  consideration  must  be  given  to  a  choice  in  displacement  sensors  since  the  CD  involves  strong 
electromagnetic  fields  during  the  discharge.  Such  fields  will  interfere  with  any  locally  mounted  electronic  force 
sensors,  such  as  piezo  based  load  cells.  This  requires  that  a  non-electronic  or  remote  sensing  sensor  be  used.  This 
interference  is  only  a  factor  during  the  thruster  firing.  Thus,  for  sensors  that  recover  from  such  interference  quickly 
the  effect  is  negligible  for  Isp  measurements.  Linear  voltage  differential  transformers  (LVDTs)  have  been 
successfully  implemented  for  previous  iterations  of  similar  thrust  stands.7,8  These  sensors  have  an  advertised 
linearity  of  <  0.25%.  For  a  2.5  mm  full  range  device,  this  translates  to  a  linear  accuracy  of  better  than  6.25  pm.  The 
linear  range  and  accuracy  of  the  detector,  for  a  given  thrust  stand  configuration,  directly  set  the  maximum  impulse 
and  minimum  mass  loss  able  to  be  resolved  by  the  thrust  stand.  Unfortunately,  LVDTs  have  an  inherent  excitation 
frequency  on  the  order  of  103  Hz,  which  limits  the  response  time  and  makes  it  unfeasible  for  time  accurate  profile 
measurements  of  the  CD,  which  has  a  pulse  width  on  the  order  of  300  ps.5  For  the  time  accurate  profile 
measurements,  one  alternative  to  the  LVDT  is  an  interferometric  system.  One  particular  laserinterferometric 
vibrometer  has  an  advertised  displacement  range  of  ±20  mm,  a  resolution  on  the  order  of  0.3  nm  and  a  sampling  rate 
of  up  to  106  Hz.’’  In  addition,  an  interferometric  system  is  not  influenced  by  the  EMI  given  that  the  working 
distance  of  the  interferometer  allows  the  electronics  to  be  mounted  outside  the  test  cell. 

In  order  to  correlate  position  measurements  with  known  forces  and  impulses,  the  stand  must  be  calibrated  in  the 
same  configuration  that  the  thruster  testing  is  conducted.  The  most  straightforward  system  is  to  place  calibrated 
weights  on  the  stand  to  create  a  known  force.  Unfortunately,  this  method  is  not  easily  replicated  for  in-situ  vacuum 
calibration.  Previous  iterations  of  thrust  stands  have  used  electrostatic  combs  to  generate  repeatable  and  known 
attractive  forces.10  Unfortunately,  the  maximum  steady  state  force  generated  by  such  a  system  has  only  been 
demonstrated  on  the  order  of  10"6-10~3  N.  With  a  thrust  stand  period  on  the  order  of  1  second,  impulses  on  the  order 
of  10"7-10"4  N  s  are  all  that  a  comparably  sized  system  could  impart  and  still  adhere  to  a  pulse  width  much  less  than 
the  stand  period.  This  is  too  small  for  an  expected  CD  impulse  magnitude  of  10’1  N  s5.  To  cover  the  lower  region  of 
impulse  calibration  and  allow  for  steady  state  force  calibration,  a  larger  scaled  electrostatic  comb  system  has  been 
designed  and  built  to  impart  up  to  5x1 0~2  N  of  force  to  the  thaist  stand.  To  cover  the  higher  region  of  impulse 
calibration,  a  piezoelectric  impact  hammer  will  be  implemented  to  calibrate  in  the  lO^-lO1  N  s  range.  Such  a  system 
offers  a  much  higher  range  of  calibration  values,  but  may  not  be  easy  to  implement  in  a  vacuum  environment. 


' '  SIOS,  SP-Series  plane-mirror  miniature  interferometer,  http://sios.de/ENGLISCH/PRODUKTE/SP_E_FOT.PDF 
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Figure  3  Experimental  thrust  stand,  calibration  system,  LVDT,  and  magnetic  damping. 

When  choosing  a  set  of  flexure  pivots  for  the  stand,  the  primary  concern  for  the  specific  impulse  measurement  is 
that  the  change  in  beginning  and  ending  positions  due  to  CD  mass  loss  is  large  enough  to  be  resolved  by  the 
displacement  sensor  and  ancillary  acquisition  hardware.  For  the  steady  state  condition  before  and  after  the  CD 
discharge,  the  displacement  of  the  thrust  stand  is  only  dependant  on  the  magnitude  and  position  of  the  mass  loss  and 
the  K  value  of  the  stand.  Figure  4(A)  shows  the  expected  LVDT  reading  on  the  designed  thrust  stand  shown  in 
Figure  3  for  a  12  mg  mass  loss,  the  amount  predicted  by  the  CD  model  in  Ref.  6.  Given  a  data  acquisition  hardware 
resolution  of  better  than  1  mV,  any  of  the  spring  constants  in  the  plotted  range  would  be  able  to  resolve  a  single  shot 
mass  loss  measurement,  with  a  preference  for  larger  displacements  going  to  smaller  spring  constants.  In  addition  to 
the  displacement  due  to  mass  loss,  the  spring  constant,  as  well  as  the  moment  of  inertia,  affects  the  maximum 
displacement  range  of  the  thrust  stand  for  a  given  impulse.  The  displacement  range  decreases  with  increasing  spring 
constant  and  moment  of  inertia.  Since  there  is  a  diminishing  return  for  adding  mass  to  the  thrust  stand,  a  preference 
for  smaller  displacement  ranges  within  the  range  of  the  displacement  sensor  gives  preference  to  higher  spring 
constants.  This  leads  to  a  trade  that  yields  an  optimum  value  for  the  spring  constant  when  balancing  between  the 
sensitivity  needed  for  the  mass  loss  measurements  and  the  restoring  force  needed  for  the  total  impulse 
measurements.  Figure  4(B)  shows  the  effect  of  the  moment  of  inertia  on  the  maximum  displacement  range  of  the 
thrust  stand  for  several  values  of  the  spring  constant. 

0.12  250 


Thrust  5tand  5pring  Constant  (N  m  /  rad) 


(A)  (B) 

Figure  4  (A)  Thrust  stand  deflection  vs.  K  and  (B)  thrust  stand  deflection  range  vs.  /  for  multiple  K . 


IV.  Results  and  Discussion 

As  a  iterative  first  step  to  Isp  measurements,  a  91.44  cm  thrust  stand,  calibrated  by  electrostatic  combs  and  sensed 
by  a  2.5  mm  full  range  LVDT,  was  built.  This  iteration  represented  only  the  elements  required  for  a  working  thrust 
stand.  The  model  coefficients  for  the  thrust  stand  can  be  seen  in  Table  1.  The  measured  K  value  matches  the 
advertised  unloaded  spring  constant  from  the  manufacturer  (published  as  ±10%)  within  1.5%.  Given  these  values 
and  the  expected  CD  impulse,  this  iteration  had  a  moment  of  inertia  which  would  yield  an  LVDT  deflection  range, 
from  Figure  4,  much  larger  than  the  useful  range  of  the  sensor.  Changing  the  magnets  or  torsion  springs  (thus  C  or 
K)  yields  a  discrete  set  of  combinations.  However,  altering  the  moment  of  inertia  offers  near  continuous  possibility 
by  simply  adding  mass  to  the  rotational  arm.  Thus  the  simplest  way  of  reducing  the  deflection  range  for  a  given 
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impulse  is  to  increase  the  mass  moment  of  inertia.  The  resultant  characteristics  for  the  optimized  stand  are  also 
given  in  Table  1. 


Thrust  Stand 

1st  Iteration 

Operational 

K 

41.75 

36.22 

c 

0.191 

0.428 

I 

0.245 

1.988 

Table  1  Thrust  stand  characteristics 

As  verification  for  the  thrust  stand  model  discussed  in  the  previous  section,  equations  (3)  and  (4)  were  solved  for 
an  experimentally  tested  7.25x1 0-4  N  s  calibration  impulse.  The  comparison  between  the  experimental  first  iteration 
and  operational  thrust  stands  and  the  model  can  be  seen  in  Figure  5,  with  the  model  shifted  to  account  for  an 
arbitrary  DC  offset  in  the  experiment.  Note  that  the  experimental  and  the  model  line  data  almost  completely 
coincide.  The  change  in  I  and  observed  change  in  K  are  properly  reflected  in  a  change  in  measured  period  and 
maximum  deflection  range  between  the  two  iterations. 


—Experiment 
— Model 


—Experiment 
— Model 
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(A)  (B) 

Figure  5  Model  comparison  between  (A)  first  iteration  and  (B)  operational  thrust  stands  for  a  3.63  mNs 

impulse. 

In  order  to  test  the  CD  attachment,  data  acquisition  and  control,  and  the  power  transfer  connections,  the  CD  was 
attached  to  the  first  iteration  thaist  stand.  Because  the  CD  at  operational  power  levels  would  result  in  an  impulse 
larger  than  the  first  iteration  thrust  stand  would  be  able  to  measure,  the  CD  was  operated  at  very  low  power  levels. 
The  displacement  range  measured  by  the  LVDT  on  this  stand  versus  the  voltage  potential  on  the  discharge  capacitor 
bank  can  be  seen  in  Figure  6.  From  Ref.  6,  the  total  impulse  should  increase  with  an  increase  in  capacitor  potential, 
a  trend  that  is  reflected  by  the  increase  in  measured  deflection  range. 

i 


500  700  900  1100  1300  1500 

Capacitor  Potential  (V) 

Figure  6  Thrust  stand  displacement  range  vs.  CD  capacitor  potential. 

While  operating  at  its  full  discharge  potential,  the  CD  is  expected  to  produce  ~10_1  N  s  of  impulse  with  a 
propellant  mass  of  12  mg6.  The  expected  output  of  the  LVDT  for  such  a  firing  on  the  operational  thrust  stand  can  be 
seen  in  Figure  1.  This  represents  the  operational  levels  for  which  the  thaist  stand  was  designed  and  optimized.  If 
experimental  values  are  found  to  be  significantly  different  than  those  predicted  by  the  CD  model,  the  thrust  stand 
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maintains  the  flexibility  to  be  easily  optimized  for  a  large  range  of  total  impulse  values  by  variations  in  the  mass 
moment  of  inertia,  spring  constant  or  damping  coefficient. 
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